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ABSTRACT 


A new gas-contvolled heat pipe design is 
desovidyed which uses a liquid matrix reservoir ^ 
or sponge, to replace the standard gas reservoir. 
Reservoir volume may be reduced by a factor of 
five to ten for certain gas-liquid combinations, 
while retaining the same level of temperature 
control. Experiments with ammonia, butane, and 
carbon dioxide control gases with methanol 
working fluid are discussed. 
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cross-sectional area of heat pipe wall, cm 

3 

= liquid phase noncondensable gas concentration, g-moles/cm 

3 

vapor phase noncondensable gas concentration, g-moles/cm 

3 

pretransient dissolved gas concentration, moles/cm 

2 - 

gas diffusivity in liquid, cm /sec 
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Ostwald coefficient 

fraction of condenser filled with wick/fluid 
fraction of reservoir filled with wick/ fluid 
surface tension, dynes/cm 
gap width, cm 

fraction of wick/fluid composite filled with 

fraction of wick/fluid composite filled with 

See equation 6 

See equation 6 

See equation 19 

diffusive time constant, sec 

thermal time constant, sec 

average value of C 


composite 

composite 


fluid in condenser 
fluid in reservoir 
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1.0 INTRODUCTION 

This report describes analysis and experiments for a new type of 
thermal control device called a gas absorption reservoir heat pipe. 

Briefly, the device maintains a relatively stable heat source temperature 
by a stagnant gas leg in the condenser which moves in response to varying 
heat input, as occurs with a gas-controlled heat pipe. The device differs 
from the standard gas controlled heat pipe in that the gas reservoir 
associated with the heat pipe is replaced by a much smaller liquid reservoir. 
For a number of gas/fluid combinations, it can be shown that it is 
volumetrically more efficient to store gas as a dissolved solute than 
dispersed as a gas in a vapor reservoir. In the experiments described 
in this report, temperature stability has been improved by a factor of 2.5 
or more over a gas reservoir of equal volume with the gases ammonia and 
butane in methanol. For the heat pipe tested, equivalent temperature 
control would have required a gas reservoir over 5 times the volume 
of the liquid sponge reservoir. 

The absorption reservoir heat pipe has some drawbacks, however, 
including potentially longer response times to step-changes in heat transfer 
rate. Section 2, Analysis, discusses absorption reservoir design in some 
detail. Section 3, Experimental Tests, describes experimental methods, 
presents experimental data, and correlates the data with models presented in 
Section 2. Section 4 summarizes the analytical and experimental program, 
and suggests future development work in this area. 
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2.0 ANALYSIS 

The first section presents basic equations describing the steady-state 
response of a gas-absorption reservoir heat pipe and compares the derived 
equations with those for a gas-controlled heat pipe embodying a gas reservoir. 
The next section discusses possible gas-liquid combinations for room- temperature 
and cryogenic operation found from a literature search. The third and last 
section addresses the effect of rapid source strength variations, which may 
temporarily drive the gas absorption reservoir heat pipe to higher than 
equilibrium operating temperatures because of a rate-limiting process of gas 
diffusion into the liquid reservoir. 

2. 1 Basic Operating Theory 

Figure 2.1 presents a comparison between the standard gas-controlled 
heat pipe and the gas absorption reservoir heat pipe. The standard design 
comprises a condenser active length Z and a gas reservoir of volume V-. 

The equivalent gas absorption reservoir heat pipe has a liquid storage volume 
storing volume of fluid and gas, but is otherwise identical to the first 
heat pipe. Each has a total cross-sectional area A^ in the condenser, which 
is independent of axial position. 

If for purposes of comparison the gas front is assumed sharp (i.e., a 
flat front), then simple closed-form models of the two systems are possible. 

For thermal comparison, assume it is necessary to design a gas-controlled heat 
pipe with the following constraints: 

1. At vapor core temperature T.j, the gas front shall be at Jig = Jl^. 

2. At vapor core temperature T 2 , the gas front shall be at Jig = 0. 

3. The sink temperature is Tg, and the gas zone can be assumed 


to also be at Tq. 


Q- 

NOTE: 



STANDARD GAS 
RESERVOIR 


GAS ABSORPTION 
RESERVOIR 


FIGURE 2.1 A COMPARISON OF GAS RESERVOIR AND GAS ABSORPTION RESERVOIR 
HEAT PIPE CONDENSER ASSEMBLIES 
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4. The gas reservoir temperature is at Tq, and the reservoir is wicked. 

5. Ambient temperature is constant. 

The wicked reservoir heat pipe described above operates nominally at 
0.5 (T^ + T^) with control range ±0.5 (T 2 - T^). By invoking conservation of 
noncondensable gas (NCG) at constraints (1) and (2), the NCG reservoir can be 
sized in terms of other thermophysical characteristics: 


Vr 


( ^Sl ~ ^so) 
(^S2 ■ ^Sl) 




( 1 ) 


where corresponds to the saturation vapor pressure at temperature T^. , and 
3^. and Bp are the fractions of condenser and reservoir volumes respectively 
filled with wick and fluid. 

If the gas reservoir is replaced by a liquid- vapor matrix with total 

volume and liquid volume UpB^V^^p, then under conditions of vapor-liquid 
equilibrium, the concentrations of NCG in the two phases will be related by 


C 


9 ^ 



2 

(g-moles/cm ) 


( 2 ) 


where is the molar gas density in the liquid phase and Cg^ is the molar gas 
density in the vapor phase. The factor a is the Ostwald coefficient and is a 
constant for dilute solutions, although it may be dependent on concentration at 
high concentration and/or high gas partial pressure. The fraction of reservoir 
volume filled with fluid is n^B^; is the fluid volume fraction of the reservoir 
wick/fluid composite. The fractional fluid fill in the condenser wick/fluid 
composite is correspondingly given by If conservation of NCG is invoked at 
conditions (1) and (2), the following is obtained: 
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where 


"o " [Wc ^ ''jA] (Condition 1) (3) 

^0 " (Condition 2) (4) 


Cgv<^> 



(5) 


and 

If expressions (3) and (4) are set equal to each other and is solved for, 
it is found that 



■ ^Sl^ 


A I 


c c 



(7) 


If, for purposes of comparison, it is assumed in equation (7) that is 
very close to 1.0 and 3^ is much less than 1.0, then it is clear that expres- 
sions (1) and (7) differ by only the multiplicative constant l/a3^ if 3^. and 
3^ are both small in equation (1) and a is much larger than one. That is, 
the required liquid reservoir volume for a given amount of temperature control 
is on the order of l/a3^ times the gaseous volume needed in a typical gas 
reservoir heat pipe. 

Phrased another way, if gas/fluid combinations can be found where a is 
significantly greater than 1.0 and is small, then it is volumetrical ly 
more efficient to store the gas as a dissolved solute than dispersed in a vapor 
phase reservoir. Such combinations are possible but common control gases do 
not satisfy this criteria. 
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2.2 Candidate Gas/Fluid Systems 

As discussed in Section 2.1, it is volumetrical ly more efficient to store 
a noncondensable gas in the form of a dissolved solute than to store the gas 
dispersed in a vapor phase if the Ostwald coefficient is greater than 1. It is 
thus desirable to use systems with the highest Ostwald coefficient possible for 
gas absorption heat pipes, unless very high solubility contributes to deleterious 
side effects. 

Stable heat pipe operation has been previously demonstrated for a in 
the range 20 to 30. In Reference 1, Saaski operated both Cl 2 /CCl^ and CI 2 /CFCI 2 
systems in the reflux mode for temperatures between 10° and 100°C with satisfactory 
component separation, Ostwald coefficients as a function of temperature are 
shown in Figure 2.2. Both carbon tetrachloride and R-11 are unfortunately 
rather poor heat pipe fluids, and a literature search was accomplished to 
determine if there existed more optimum fluid/gas combinations. 

Table 2.1 presents Ostwald coefficient data from Reference 2 for a large 
number of gas/liquid combinations at 25°C. The common control gases such as He, 

N 2 and Ar are seen to have Ostwald coefficients strictly less than 1, and are 
unsatisfactory. Large molecules such as Kr, CO 2 , Xe and C 2 Hg do have Ostwald 
coefficients greater than 1, but the largest value is on the order of 5. This 
is somewhat marginal. To demonstrate a significant advantage over present gas 
reservoir designs, it would be desirable to use a gas/fluid combination with 
an Ostwald coefficient of 10 to 20. As discussed in Section 2.3, some void 
volume must be left in the reservoir for optimum reservoir response, and a 
fluid/wick volume fraction of 0.5 is a realistic goal. Reservoir size 
reduction would therefore be on the order of 1/5 to 1/10. 


OSTWALD COEFFICIENT 
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Y (10‘^K) 


FIGURE 2.2 THE SOLUBILITY OF CHLORINE IN CCl^ AND R-11 (CFCI 3 ) 



TABLE 2.1. EXPERIMENTAL OSTWALD COEFFICIENTS AT 25°C FOR A NUMBER OF GASES AND FLUIDS 
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Table 2.2 presents Ostwald coefficients for several room temperature 
combinations. It can be seen that several gas/liquid combinations are available 
for laboratory ambient heat pipes, and system solubilities are in a useful 
range. In some cases binary combinations such as water-methanol could be more 
appropriately classified as two-fluid media because of an extremely high mutual 
solutility. In these situations, the Ostwald coefficient is not quantitatively 
useful because the solubility is dependent on both pressure and temperature, 
and dissolution of large amounts of solute leads to nonideal behavior. 

Methanol was selected as a working fluid for room temperature tests. 

In addition to being widely used as a heat pipe working fluid, a relatively 
low vapor pressure minimizes gravity-stratification of the condenser gas front. 
Non-perpendicular gas fronts have been observed in high-pressure ammonia- 
xenon heat pipes during 1-g tests, and this behavior has been observed to create 
erratic response. By using a low pressure (corresponding low density) 
working fluid having relatively low molecular weight, stratification with 
control gases of dissimilar molecular weight is minimized. Figure 2.3 shows 
experimental data from Reference 5 on the solution of various gases in methanol 
as a function of temperature. In all cases, the gas partial pressure is 1.0 
atmosphere. In experiments discussed in Section 3, the gases nitrogen, carbon 
dioxide, butane, and ammonia were used. Ammonia is not presented in Figure 2.3, 
but the solubility is estimated to be higher than SO 2 . In all cases, the 
solubility is approximately exponential in temperature. 

This exponential behavior is also found for highly soluble gases in 
cryogenic working fluids. Figures 2.4 and 2.5 show data for several cryogenic 
combinations. References for these combinations are given in Table 2.3. 
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TABLE 2 


Solvent 

Hexane 

Benzene 

Benzene 

Menthanol 

Menthanol 

Menthanol 

Menthanol 

Menthanol 

Water 

Water 

Water 


.2 ROOM TEMPERATURE GAS-LIQUID COMBINATIONS HAVING 
HIGH SOLUBILITY 


Solute 

Temperature 

(°c) 

Ostwald 

Coefficient 

Ref 

n-Propane 

25 

23.6 

(3) 

n-Propane 

25 

16.0 

(3) 

n-Pentane 

16 

312.0 

(4) 

Propane 

25 

3.4 

(5) 

Carbon Dioxide 

12 

4.1 

(5) 

Butane 

12 

28.0 

(5) 

Sulfur Dioxide 

25 

83.0 

(5) 

Carbon Dioxide 

59 

39.0 

(6) 

Ammonia 

25 

40.7 

(5) 

Sulfur Dioxide 

25 

34.0 

(5) 

Menthanol 

100 

254.0 

(7) 


SOLUBILITY, MOLE FRACTION 
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2.3 2.5 2.7 

log^q(t) (t in K) 

FIGURE 2.3 SOLUBILITY OF VARIOUS GASES IN METHANOL AS A FUNCTION OF 
TEMPERATURE (FROM REFERENCE 5) 










(OSTWALD COEFFICIENT) 
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TABLE 2.3. REFERENCES FOR CRYOGENIC BINARY 
GAS/LIQUID SOLUBILITY DATA 


Gas/Fluid System 
He/H2 
He/Ar 
He/N2 
He/02 

H2/CH4 

Ne/Ar 

Ne/N2 

N2/O2 

N2/CH^ 

Ar/CH^ 

CH4/C2Hg 


Reference Numbers 
8 

9, 10 
11 
10 
12 
9 
11 

13, 14 
15, 16 
15, 16 
7 
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Table 2.4 gives fitting constants A and B for cryogenic combinations 
of interest where the Ostwald coefficient is modeled as 

« ^ A EXP(+B/T) (8) 

and where T is the system temperature in degrees Kelvin. The exponential 
coefficient B ranges from 550 to 1850 K and the pre-exponential factor is 
between 0.07 and 0.10. Systems of potential interest include 

1 . methane/ethane 

2. argon/methane 

3. nitrogen/methane 

4. nitrogen/oxygen 

The solubility of helium in hydrogen is surprisingly low as shown in Figure 2.5 
and, hence, this technique would not be practical for liquid hydrogen nor, of 
course, helium. 

The general exponential dependence of solubility on temperature is of 
both beneficial and adverse effect. Because of the rapid solubility decrease 
with temperature, condensate originating near the cold gas front becomes 
supersaturated as the fluid warms in transit to the evaporator, and gas stripping 
is augmented. This serves to maintain a sharp gas front. Conversely, the 
sensitivity to temperature can create problems in set-point maintenance. If 
liquid reservoir temperature is not maintained constant, the amount of dissolved 
gas in the liquid will fluctuate in response to temperature variations. The 
most direct consequence of this will be a fluctuation in both gas front position 
and heat pipe operating temperature. This sensitivity to reservoir temperature 
is also found in the wicked reservoir heat pipe where the reservoir working 
fluid partial pressure is exponentially dependent on reservoir temperature. 
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TABLE 2.4. CURVE-FITTING CONSTANTS FOR BINARY CRYOGENIC 

SYSTEMS SHOWING EXPONENTIAL SOLUBILITY BEHAVIOR 


Gas/Liquid System 

A 

B (K) 

CH 4 /C 2 Hg 

0.0713 

848 

Ar/CH^ 

0.0748 

673 

N 2 /CH 4 

0.0821 

549 

N 2 /O 2 

0.100 

599 


a 
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2.3 Transient Response-Reservoir Design 

Two criteria dominate the physical design of an absorption reservoir. 
First, there must be adequate contact between the absorbing liquid and the 
vapor phase to ensure speedy equilibration of dissolved gas content in response 
to changes in thermal input because of typically low solute diffusion rates in 
liquid. Second, the vapor spaces used for vapor communication cannot be so 
small that they prime and drain liquid away from working sections of the heat 
pipe. In addition, complete filling of the reservoir with liquid will result 
in slow transient response. 

In a representative reservoir design, the liquid would be contained in 
a very porous sheet of capillary material wrapped in a spiral so that the 
principal fluid priming radii of curvature for the inter-spiral gas gaps would 
sum to a constant. 



where r.j and r^ are principal radii of curvature for the gap, and y is surface 
tension. 

For a spiral of constant pitch and inter-spiral gap width 6, this 
pumping pressure is approximately 


AP 


6 


2l 

6 


(10) 


If it is necessary to prime a capillary structure such as an artery of 
radius r^, then to ensure proper priming of the artery it must be true that 

a 

> lx 
r, 6 


( 11 ) 
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If the heat pipe is tested in 1-g then the effects of gravity must also be 
included. 

Equation 11 presents a criterion relating minimum size of reservoir gas 
gaps to characteristic capillary pumping structure dimensions in the heat pipe. 
To characterize maximum fluid absorber thickness, and thereby define the void/ 
volume ratio in the reservoir, it is necessary to establish reservoir response 
to fluctuating heat inputs. 


2.3.1 Reservoir Response-Slab Model 

A representative concentration profile for a section of absorber is 
shown in Figure 2.6 after a step-decrease in power, i.e., a reduction in non- 
condensable gas partial pressure. If a unit cell of the absorber is modeled 
as a thin disc or slab, then the average time-dependent excess concentration 
is given as the series 


— * / X 

C (t) 


_8_^ EXP ( - (2N+l)Vt) 
tt^N=0 (2N+1)^ 


( 12 ) 


where D* = (J) 

and Dgjj is the binary gas/fluid diffusivity and I is the overall thickness of 
a given disc. A derivation of this equation is given in Appendix A. The excess 


concentration C is defined as 

^ ^ <^ - ^o> 


(13) 


where C^is the pre-transient dissolved gas concentration and is the post- 
transient equilibrium gas concentration. This dimensionless concentration is 

2 

shown as a function of the grouping D^j^t/S. in Figure 2.7. For this grouping 
greater than 0.025, the dimensionless concentration is given approximately as 
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♦-X 




FIGURE 2.6 ONE-DIMENSIONAL SLAB MODEL OF ABSORPTION RESERVIOR UNIT 
CELL. WITH A STEP-CHANGE IN EXTERNAL GAS CONCENTRATION, 
THE POROUS LIQUID-SATURATED ABSORBER RE-ESTABLISHES 
EQUILIBRIUM BY GAS/LIQUID DIFFUSION. 



DIMENSIONLESS CONCENTRATION (c 
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C* s 0.82 EXP (-9.89 Dt/£^) (14) 

★ 

At time equal zero, C =1., while at t = «>, c = 0. At some finite 
time C is negligibly small. If c = 0.01 is arbitrarily selected as negligibly 
small, then a characteristic time constant t for diffusion-driven equilibrium 
is obtained 

T = 0.446 (secs) (15) 

Because cryogenic heat pipe design is currently of high interest and 
fluid diffusivities can be expected to be low, values of t have been calculated 
for several cryogenic fluids using the following expression for the diffusivity 

9. 76(10"^) 

“ab = TS^ 06) 

y. a 

b a 

where 

2 

Dgjj = diffusion coefficient, a->-b (cm /sec) 

= solvent molecular weight 

ip = association parameter (water = 2.6, methanol = 1.9, benzene = 1.0, etc. 

T = temperature (°K) 

Ujj = solvent viscosity (centipoise) 

= molecular diameter, A° 

a 

This is a modified form of the Wilke-Chang empirical diffusivity formula 
as given in Reference 2. Since this was derived from room temperature data, it 
was necessary to check the model against experimental cryogenic data. This was 
done for the data of Ricci (Ref. 17) for various gases in liquid nitrogen. 

Figure 2.8 shows correlation of Ricci's data with the empirical model (16) for 
the gases Ne, Ar, and CH^. The temperature dependence of diffusivity is 
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represented very well, although the absolute values differ by as much as 30%. 
However, this is adequate agreement for scoping purposes. 

In the previous section, several cryogenic gas/liquid combinations were 
identified as of potential interest. These were 

1 . methane/ethane 

2. argon/methane 

3. nitrogen/methane 

4. nitrogen/oxygen 

With the model (16), the diffusivities shown in Figure 2.9 were calculated. 

-52 -52 

The diffusivities are typically from 2 x 10 cm /sec to 10 x 10 cm /sec. 

Before using these values of 0^^^ to calculate t, the values were multiplied by 
0.40 to account for permeability of the porous slab used as a gas absorber. 

This correction factor was based on the experimental data of Saaski (Ref. 2) 
for double- thickness square-weave screens. 

Characteristic time constants are presented in Table 2.5 for the fluids 
previously noted, using Equation 16 and a correction factor of 0.40 to 
A slab thickness of 0.050 cm was selected as representative of absorber thickness 
the time constant scales as Z for other thicknesses. For this particular 
choice of £, time constants of 0.6 to 2.3 minutes are indicated. The time 
constants are therefore not excessively long, but care must be exercised in 
design so that the characteristic fluid films are of the order of 0.050 cm in 
thickness. Even an increase to only 1.0 mm thickness results in time constants 
of 2.4 and 9.2 minutes, respecti vely; the latter equilibration time could be 
unsatisfactory in many thermal control systems. 


_oi 
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4 6 8 10 12 14 16 18 


Y O0‘^/K) 

FIGURE 2.9 PREDICTED DIFFUSIVITIES FOR VARIOUS GAS-LIQUID COMBINATIONS 
THAT ARE OF INTEREST FOR GAS-ABSORPTION RESERVOIRS AT 
CRYOGENIC TEMPERATURES 
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TABLE 2.5. 

DIFFUSION TIME CONSTANTS FOR VARIOUS 
LIQUID BINARY SYSTEMS U = 0.050 cm) 

CRYOGENIC GAS/ 


Temperature 

Time Constant 

Gas/Liquid 

(°K) 

(minutes) 


86. 

2.3 

N^/CH^ 

100. 

2.0 


140. 

0.7 

Ar/CH^ 

100. 

1.6 


140. 

0.6 

CH^/C^Hg 

148. 

1.8 


213. 

0.6 
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2.3.2 Coupled Heat Transfer-Diffusion Problem 

Section 2.3.1 gives an order-of-magnitude model for diffusion times 
associated with absorption reservoir heat pipes. More precise modeling of 
an actual device requires coupling of heat transfer processes in the heat 
pipe with diffusion processes in the reservoir. This section describes 
analytical and numerical modeling of such a system. 

2. 3. 2.1 Conservation of Noncondensable Gas 

Unsteady-state modeling of heat pipe response requires coupled 
equations for conservation of energy and conservation of noncondensable 
gas. At time t < 0 the heat pipe is in equilibrium at a power level 
and the total gas content is given by 

% = [Wg"”" Vr] (17) 

where C (0 ) denotes the gas concentration for t < 0, and Si (0)is the 

y q 

initial length of the gas leg. In this representation, the gas leg is 
assumed to end abruptly at the transition to free-flowing vapor, i.e., a 
'flat front' is assumed. At t = 0, a step change in power level to 
occurs, and conservation of noncondensable gas requires 



N = 
0 

c (t) Fa 0 £ (t) + v„ e’*’! 

gv' 'I c c g' ' £r r 

(18) 

where 

< = 

1 + (otn^x^-1) 

(19) 

and 

X = 

r 


(20) 
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It is assumed that diffusion of dissolved gas into and out of the 
condenser wick is very rapid, and that the primary non-equilibrium condition 
occurs in the reservoir fluid. Deviation from equilibrium conditions at 
t > 0 is represented by x^, which is the ratio of actual dissolved gas 
concentration to the dissolved concentration at the vapor-liquid interface. 
This interface is always assumed to be in equilibrium, that is. 


Cgj^(t) = aCg^(t) (interface condition) (21) 

The average dissolved gas concentration in the reservoir is 
<Cgjj^(t)>. For a slab reservoir of thickness Z, the average dissolved 
concentration is expressed as 





( 22 ) 


This assumes that the slab is in vapor contact on both faces, as in Figure 2.6. 


2. 3. 2. 2 Conservation of Energy 

At any instant in time, heat transferred to the heat pipe must be 
accounted for by heat transfer through the condenser and/or by a temperature 
increase of system thermal mass. For modeling ease, system thermal mass is 
assumed always to be in equilibrium with the vapor core temperature, T^(t). 
Conservation of energy at t > 0 results in 







(Vs) 


mC 


( 23 ) 
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Heat loss via the condenser is the sum of a conductance loss from 
the unblocked condenser zone a conductance loss from the 

blocked-off zone, In the latter case, the heat pipe wall in the blocked- 

off zone acts as an extended fin coupled to the sink with a per-unit-length 
conductance G^. At the fin root, i.e., at the gas-vapor flat-front inter- 
face, the wall is assumed to be uniformly at temperature T^. For a long 
condenser, where the gas zone simulates an infinitely long fin, the 
equivalent fin-conductance is 

G = <24) 


2. 3. 2. 3 Numerical Methods 

In the previous sections, the necessary conservation laws have been 
outlined: At any point in time, both noncondensable gas and heat input 

must be accounted for. However, the equations are coupled through the gas- 
leg length, which in turn is dependent on the amount of dissolved gas 
in the liquid reservoir. 

The various equations (18) through (23) have been solved simultaneously 
by numerically solving the diffusion equation in the slab reservoir wick 
and step-wise integrating equations (22) and (23). At each new time-step, 
was iterated until T^, £^, and dC^^/dx did not change by more than an 
arbitrary e per iteration. When a consistent set of variables occurred, 
the integration in proceeded to the next time-step. The trapezoidal rule 
was used for integrating both dC^^/dx and dT^/dt. The iterative sequence at 
each time step was 
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1. Calculate from equation (23) using previous value of £ ; 
calculate 

2. Using and equation (21), calculate dC^j^/dxi using 

finite differences. ^ ^ 

3. Calculate the average dissolved gas concentration using (22). 

4. Determine 0^^. from equations (19) and (20). 

5. Determine a new from equation (18). 

6. Compare with 2,^^. of previous iteration. If A£g < e, 
continue to next time step. If A£g > e, go to step 1 and 
perform calculations using updated 2,^^ . 

Representative unsteady-state solutions for the time-dependent vapor 
core temperature T^(t) are given in Figure 3.7. These solutions are based 
on the experimental device detailed in Section 3.1, Experimental Methods 
and Apparatus. 

The numerical solutions are useful in quantitatively predicting 
performance for a well -documented heat pipe system. However, they are also 
useful in establishing the relative importance of thermal time constants and 
diffusional time constants. That is, from equation (23) a thermal time 
constant can be derived 


mCp 

■^th " Gp<£p-£g(t)> 


(25) 


which is based on an average nonblocked condenser length <£p-2g(t)>. The 
ratio of thermal to diffusion time constant is 


^ = '"^p^ab 

^D Gp£^<£p-2,g(t)> 


( 26 ) 
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The numerical solutions indicate that an overshoot in the curve 
T^(t) occurs for x'*' on the order of 20. For x'*’ less than 20, the heat 
pipe temperature temporarily rises above the final steady-state operating 
temperature because of the lag in gas permeation of reservoir wicking. This 

criterion for x^ should be a valid rule-of-thumb for estimating absorption 
reservoir response for diverse heat pipe designs and gas/fluid selections. 
However, actual experimental data do not indicate an effect of diffusion 

limiting absorption for x^ much less than 20, and this divergence in 
behavior remains to be explored. 

2.4 Analytical Summary 

The temperature-control capability of a gas absorption reservoir heat 
pipe has been derived for both steady-state and transient conditions. 

The control band of a given heat pipe design may be reduced by up to 
a factor of 10 by replacing a gas reservoir with an equal-volume absorption 
reservoir. However, absorption reservoir design requires careful consideration 
of mass transfer, as the dissolved gas content within the reservoir should 
respond rapidly to changes in heat input to reduce temperature fluctuations 
associated with changes in power level. Numerical models of unsteady-state 
behavior indicate that the ratio of heat pipe thermal time constant to 
diffusional time constant must be greater than or equal to approximately 20. 
to minimize temperature over-shoot or under-shoot when a step-change in 
power occurs. The temporary rise or fall in temperature of an absorption 
reservoir heat pipe can occur because gas diffusion in the reservoir liquid 
does not keep pace with the rate of temperature change based on system 
thermal mass. 

A literature search has resulted in a number of potential gas/liquid 
combinations suitable for gas absorption reservoirs. Combinations have been 
found for both room temperature and cryogenic temperature ranges. 
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3.0 EXPERIMENTAL VERIFICATION 

Experimental test methods are described in Section 3.1. In Section 3.2, 
absorption reservoir experimental results are presented, while data are cor- 
related using models described in Section 2 in Data Interpretation, Section 3.3. 

3. 1 Experimental Methods and Apparatus 

3.1.1 Experimental Design 

Two complementary test systems have been used for experimental verifi- 
cation of absorption reservoir performance. Both systems use the same basic 
heat pipe assembly composed of evaporator, condenser, and adiabatic sections. 

The evaporator and condenser were each 15.24 cm long, with an adiabatic 
section 15.0 cm long. All three sections were 1.288 cm inside diameter. 

Nominal wall thickness in the condenser was 0.0650 cm, while the wall was 0.080 
cm thick in the adiabatic and evaporator sections. A circumferential wall wick 
of 1 layer of 200 mesh square weave 304 stainless screen was used in all 
three sections. No bypass wick was used, and the heat pipe was operated in 
a slight reflux. These design and testing simplifications ensured accurate 
estimation of fluid volume factors without inclusion of generally uncertain 
correction factors for fillets. Stainless steel was used throughout. The 
condenser heat exchanger was anodized aluminum and was isolated from the 
condenser proper by an air gap of 0.022 cm for a per-unit-length gap conduc- 
tance of approximately 0.0443 W/cm-K. Four 6-32 nylon adjusting screws 
were used at each end of the condenser block to ensure a uniform air gap. 

A slot in the wall of the condenser cooling block enabled direct spot-welding 
of the thermocouples to the condenser wall. 
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Two reservoir assemblies were used in testing. The first reservoir, 
coupled with a glass-walled section for visual observation, was formed from 
a spiral of 100 mesh stainless square-weave screen. This method of reservoir 
manufacture was found to be difficult with existing tooling and to have a 
low volumetric efficiency; that is, the fraction of reservoir volume filled 
with liquid was quite low. The glass-walled reservoir was primarily used 
for visual studies of fluid absorption; the experimental results presented 
in Section 3.2 deal exclusively with the all-metal reservoir system described 
below. 

The metal -walled reservoir used for quantitative measurements of 
reservoir function is shown in Figure 3.1. In detail (A) of Figure 3.1, an axial 
cross section of the reservoir assembly is shown. The fluid absorption 
wi eking was a single layer of 0.0635 cm thick felt-metal on the reservoir 
wall and a central cylindrical wick of 0.673 cm O.D. and 0.318 cm I.D. The 
central wick assembly is shown in dotted profile in (A) and in cross section 
in (B). The inside and outside diameters of the central wick are a single 
layer of 400 mesh stainless screen. The internal screened tubes, used to 
give the structure rigidity as well as a small priming radius, are also of 
400 mesh screen and approximately 0.062 cm in diameter. Thirteen tubes were . 
used on the inner layer and 21 tubes on the outer layer. The tubes were 
sealed at the end facing the condenser and maintained open at the other end 
to minimize bubble entrapment. 

The entire reservoir is 6.81 cm long and 1.29 cm in diameter, and has 
a 41% volumetric fluid load. Reservoir parameters needed for modeling are 
given in Table 3. 1. 
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COOLANT 

BOLT-HOLE 
0.206 cm I.D. 


RESERVOIR 
GUARD RING 
(BRASS) 



CONDENSER 
HEAT SINK 



1 

I 


&SB 

B5 




1.29 cm 


6.81cm 



METAL WICK 
(0.0635 cm) 

FLANGES 


(A.) CROSS-SECTION OF RESERVOIR ASSEMBLY 



(B.) END-VIEW OF ABSORPTION RESERVOIR CENTRAL WICK 


FIGURE 3.1 DETAILS OF ABSORPTION RESERVOIR DESIGN. WICK ON WALL WAS ONE LAYER OF 
FELT-METAL. CENTRAL RESERVOIR WICK SHOWN AS OOHED PROFILE IN (A) IS 
SHOWN IN END-VIEW IN (B). WITH ALL WICKIN6 PRIMED, 41* OF RESERVOIR 
VOLUME WAS FLUID-FILLED. 
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TABLE 3.1. RESERVOIR SPECIFICATIONS/CONDUCTANCE PARAMETERS 


Reservoir total volume 

8.87 cm^ 

Wick/fluid volume fraction (3^) 

0.465, 0.0375 

Fluid volume fraction in wick (n^)» 

0.886, 0.68 

Wick thickness 

0.068 cm 

Reservoir and heat sink temperature^ 

11 - 13°C 


0.0443 W/cmK 

Heat pipe orientation 

0.24 cm favor- 
able tilt 


^Reservoir guard ring and condenser sink both cooled by same 
water circuit. 
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3.1.2 Testing Methods 

All tests were made using methanol working fluid. The heat pipe was 
initially loaded with methanol and set to 0.25 cm reflux tilt to ensure an 
absence of fluid pooling in the condenser. Operation without noncondensable 
gas for several hours was used to determine leak tightness. The nonconden- 
sable gas was then inserted in sufficient quantity to allow gas control 
operation over the highest temperature range. Several cycles of the heat 
pipe from 10 to 25 watts appeared to stabilize heat pipe response, and this 
may be related to priming of the liquid reservoir. 

After stable response was obtained, data were gathered and gas sub- 
sequently removed for heat pipe testing in a lower temperature range. 

The absorption reservoir was used with the gases carbon dioxide, 
ammonia, and butane. After the absorption reservoir test sequence was com- 
pleted, the wi eking was removed from the reservoir and testing with nitrogen 
was used to define an experimental reference. Nitrogen has negligible 
solubility in methanol, in the context of this study. 

An operating run was also made with all gas removed and an isothermal 
heat pipe to determine heat loss through the insulation package. Heat loss 
is linear and given as 

= 0.0636 (Ty-Tg) (watts) (27) 

where T„ is vapor temperature and T^ is ambient temperature. Heat loss via 

V a 

the nylon screws to the heat sink is 

Q ^2 = 0.0226 (VTg) (watts) (28) 

Data presented in Section 3.2 have been corrected for these two losses. 


36 


3. 2 Experimental Results 

Using the testing techniques given in Section 3,1.2, Figures 3.2 - 
3.4 present the vapor-core heat-sink temperature difference, as a function 

of power level for the control gases N 2 . CO 2 , and NH^. As mentioned 

earlier, the nitrogen data were taken with the liquid reservoir wi eking 
removed and hence represent a reference gas control design of equal reservoir 
volume. The temperature control band with ammonia and butane as control 
gases is less by about a factor of 2.5 relative to N 2 for a nominal control 
range of 5 to 20 watts. The experimentally determined control bands are 
given in Table 3,2 for the gases used. The Ostwald solubility coefficient 
is also given at the nominal reservoir temperature of 12.0°C. 

3. 3 Data Interpretation 
3.3.1 Temperature Control 

Table 3.2 presents a comparison of experimental and predicted 
performance for both the absorption reservoir tests with CO 2 J and NH^ 

and gas reservoir tests with nitrogen. In general, the gases and NH^ 

had, for a given heat transfer band-width, a temperature band-width 1/2 to 
1/3 that seen with the standard gas reservoir and nitrogen. This difference 
amounts to a saving in reservoir volume of a factor of 5 to 7.5, 

The data comparison in Table 3.2 is based on AT/AQ. The heat pipe 
vapor core temperature was found to be essentially linearly dependent on Q, 
the heat transfer rate, provided that the gas front resided in the condenser. 
Therefore, the control capability has been characterized on the basis of the 
slope of the function T^(Q). In each test, the slope AT/AQ was calculated by 
a least-squares fit to the data over the linear region. The data have also 



FIGURE 3.2 TEMPERATURE DIFFERENCE BETWEEN HEAT PIPE VAPOR CORE AND CONDENSER HEAT SINK FOR THE 
CONTROL GASES n-C.H,« AND CO^. DATA TAKEN WITH ABSORPTION RESERVOIR. 
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TEMPERATURE DIFFERENCE BETWEEN HEAT PIPE VAPOR CORE ANp CONDENSER HEAT SINK FOR 
AMMONIA CONTROL GAS. DATA TAKEN WITH ABSORPTION RESERVOIR 
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TEMPERATURE DIFFERENCE BETWEEN HEAT PIPE VAPOR CORE AND CONDENSER HEAT SINK 
FOR NITROGEN CONTROL GAS. THIS DATA HAS BEEN TAKEN WITH THE ABSORPTION 
WICKING REMOVED; THE RESERVOIR IS FUNCTIONING AS A STANDARD GAS RESERVOIR. 




TABLE 3.2. EXPERIMENTAL DATA SUMMARY FOR VARIOUS GAS CONTROL TESTS WITH METHANOL 
WORKING FLUID AND A RESERVOIR TEMPERATURE OF 11 to 130C 
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been organized by the nominal vapor-sink operating temperature difference 
in the range 7.5 to 10.0 watts. The resulting data ranges, and the lengths 
of the linear regions are approximately as follows 

Region 1: O 8.75 w = 22.5°C; linear span = 2.0 to 12.5 W 

Region 2: @ 8.75 w = 37.5°C; linear span = 2.5 to 15.0 W 

Region 3: AT^^ @ 8.75 w = 50.0°C; linear span = 4.0 to 17.5 W 

The experimental values of AT/AQ in Table 3.2 represent each temperature 

region investigated. For comparison with calculated values of AT/AQ, the 
individual values have been averaged. Both the averaged and individual 
experimental values are presented in Figure 3.5 for relative comparison. 

Control capability of the absorption reservoir is in qualitative 
agreement with estimated values of AT/AQ using a flat-front model and the 
dissolution models discussed in Section 2. The experimental control ratios 
compare with theory as follows 


AT/AQ for N 2 ; 

CO^; 

C4H10; 

NH3; 


20% lower than theory 
24% lower than theory 
14% higher than theory 
100% higher than theory 


The general experimental trend is a reduced effect of Ostwald coefficient 
although performance improvements are still apparent as a increases. 

Figure 3.6 shows the experimental data compared to predicted control ratios 
for an Ostwald coefficient range from 1.0 to 1000. 

The reduced influence of ex beyond about 30 may be attributable to a 
spreading out of the gas front because of high solubility. Condenser thermo- 
couples indicated a considerably longer gas front with the three soluble 
gases than with nitrogen. However, because of differences in molecular 
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FIGURE 3.5 A GRAPHICAL ILLUSTRATION OF THE EXPERIMENTALLY DETERMINED EFFECT OF CONTROL GAS SPECIES 
ON CONTROL RATIO. THE THREE GASES, C02, C4H10. AND NH3 WERE USED WITH THE RESERVOIR 41% 
FILLED WITH FLUID. THE NITROGEN TEST WAS WITH ALL WICKING REMOVED FROM THE RESERVOIR, 
AND IT HENCE SERVES AS A REFERENCE GAS RESERVOIR DESIGN. TEMPERATURE VARIATION WITH 
HEAT INPUT FLUCTUATION WAS REDUCED BY A FACTOR OF ABOUT 2.5. 



5.0 



OSTWALD COEFFICIENT 


FIGURE 3.6 RATE-OF-CHANGE OF VAPOR CORE TEMPERATURE WITH RESPECT TO 
HEAT INPUT FOR THE THREE SOLUBLE CONTROL GASES USED WITH 
METHANOL WORKING FLUID. THE THEORETICAL INFLUENCE OF 
SOLUBILITY ON TEMPERATURE CONTROL IS ALSO SHOWN, ASSUMING 
A SHARP GAS-VAPOR INTERFACE. THE BAND STRUCTURE REFLECTS 
MAX. OPERATING TEMPERATURE RANGE. 
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weight between the gases and working fluid vapor, vertical gravity-force 
stratification of the gas front could produce the same broadening. Another 
possibility is that small amounts of insoluble gas either leak into the 
heat pipe via 0-rings or hydrogen is generated from water contained in the 
methanol. The direct effect of this will be a reduction in control ratio 
which becomes increasingly important as a increases in size. Further 
investigation of this effect is warranted. 

In summary, the temperature change resulting from a given heat input 
change has been reduced by a factor of two to three for an absorption reservoir 
of identical volume to a gas reservoir. This has been demonstrated 
with a common gas-control heat pipe working fluid, methanol. Equivalent 
performance with a wetted wall gas reservoir would require a volumetric 
increase to about 5x to 7.5x the volume of the absorption reservoir. 

3.3.2 Transient Response 

Figure 3.7 presents experimental data on time response of vapor core 
temperature to a step-increase in heat transfer rate from approximately 5 
watts to 20 watts at a nominal vapor temperature of 55°C above sink temperature. 
The vapor core temperature is nondimensional ized by defining 


AT T (t) - T 

V v' VO 


AT. 


vm 


^vm ~ ^vo 


(29) 


where T^^ is the initial vapor temperature and T^^ is the final vapor tempera- 
ture. The experimental data do not show any significant effect of reservoir 
design or control gas on transient response. However, various amounts of 
over-shoot are analytically predicted for the absorption reservoir system 
depending on the diffusion coefficient of the gas in the working fluid. The 


D=10 cm /sec 
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FIGURE 3.7 TRANSIENT RESPONSE OF THE ABSORPTION RESERVOIR HEAT PIPE TO A STEP-INCREASE IN HEAT 
INPUT. RESPONSE OF THE REFERENCE GAS RESERVOIR (i.e., N 2 GAS WITH THE WICK REMOVED 
FROM 5 W to 20 W ABSORPTION RESERVOIR) IS ALSO SHOWN FOR COMPARISON. ANALYTICAL 
MODELS FOR TRANSIENT RESPONSE PREDICT MUCH LARGER TRANSIENTS THAN WERE ACTUALLY 
MEASURED, INDICATING VERY LARGE EFFECTIVE GAS/LIQUID DIFFUSIVITIES. 
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-5 

diffusivity of ammonia gas in methanol is estimated to be 5 to 10 (10 ) 

2 

cm /sec. For this range of diffusivity, about a 40% over-shoot in dimension- 
less temperature is predicted, while only a few percent was observed. 

Butane gas also displayed minimal over-shoot and under-shoot. 

Alternate explanations for this divergence in behavior remain to be 
explored. It is possible that the high solubility of these gases in methanol 
creates a turbulent mass transfer phenomena that is much more effective than 
simple diffusion. However, if this is correct, then it must be determined 
if these turbulent transfer processes will function in zero-g, since they 
could require a gravity field. 

An alternate explanation is that the diffuse gas front found for soluble 
gas/liquid combinations is a function of both radial and axial position, so 
that the absorption process and overall heat transfer coefficient are not 
adequately described by a flat-front gas control model. Vertical stratifica- 
tion of the gas front could also create experimental deviation from analytical 
predictions . 

Viability of the absorption reservoir technique has been demonstrated 
on the basis of heat transfer. However, questions raised related to transient 
response and gravity effects should be explored in further detail to fully 
substantiate performance of the reservoir in a zero-g environment. 
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4.0 CONCLUSIONS 

Significant improvements in thermal control have been experimentally 
demonstrated for the absorption reservoir. With the same volume as a 
reference gas reservoir, temperature control has been improved by a factor 
of two to three for butane and ammonia control gases with methanol working 

o 

fluid. A smaller improvement was shown for carbon dioxide, a gas with lower 
solubility in methanol. Equivalent thermal performance from a gas reservoir 
would require a volume 5 to 7.5 times that used with the absorption reservoir 
and ammonia or butane. Detailed comparisons of absorption reservoir tests, 
standard reservoir tests with equal volume, and theoretical predictions are 
given in Table 3.2. Transient response was relatively unaffected by reservoir 
design, and heat pipe response to a step-change in heat transfer rate was 
very rapid. 

However, this responsiveness is very interesting in that the experi- 
mental heat pipe was expected to exhibit a relatively large transient 
temperature over-shoot on the basis of diffusion rate constants for the 
liquid absorption reservoir. Further analytical and experimental work is 
needed to determine the basis for this unexpected performance improvement, 
since the apparent improvement in mass transfer may be related to gravity- 
assisted convection effects not present in zero-g. 

In addition, temperature measurements on the outside of the condenser 
indicated a very diffuse gas front in the case of the soluble gases ammonia, 
butane, and carbon dioxide. The gas front was relatively sharp in the case 
of nitrogen, which is relatively insoluble in methanol and very close in 
molecular weight. A diffuse gas front and/or a front that is vertically 
stratified because of gravity effects may partially account for the improved 
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gas control, independent of the soluble gas reservoir effect; the thermal 
effects of nonperpendicular gas fronts are not well understood. It would be 
of interest to redesign the test apparatus so that the heat pipe could be 
operated in a vertical reflux position with a lower molecular weight gas 
such as ammonia. In a vertical position, the ammonia-methanol vapor phase 
mixture would effectively float on the lower density methanol vapor below. 

By carefully comparing operating characteristics in both horizontal 
and vertical modes, the effects of gravity on gas front shape could be 
determined, and transient response behavior linked to either a gas-front 
related-mechanism or reservoir mass transfer mechanism. 
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APPENDIX A 

Diffusion From a Double-sided Slab 

Given the physical system shown in Figure 2.6, it is desired to solve the 
diffusion equation in one dimension 


8C _ „ 8^C 

3t - abj^2 

A dimensionless distance and concentration are defined as 


( 1 ) 


X* = TTX/jl 

C* = (C-C )/(C-C ) 

0 'mo 

Boundary conditions are assumed to be 

@ t=0, C*=l for all X* 

(3 t=0^, C* =0 for x*=0 and x*=tt 

* 

for t-^oo, C =0 for all x* 

The diffusion equation now becomes 


3C* 

9t 


D* 


a^c* 



where 



This has the Fourier sine series solution 


C* = ^ bM EXP(-NVt) sin(Nx*) 

N=1 


( 2 ) 

(3) 


(4) 


( 5 ) 


A-2 


where 

f /* C*(t=0) sin(Nx*)dx* (6) 

^ •'o 

Using the initial concentration distribution, the solution (5) is alternatively 
given as 


* 

C = 



EXP(- (2N-l)^D*t) sin (2N-l)x* 
2N-1 


(7) 


If expression (7) is integrated in x* and divided by the interval length, then 
the average concentration is found to be 


IT 


I 

N=0 


EXP r-(2N+l)^D*t1 
(2N+1)^ 


C 


( 8 ) 
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